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Abstract 

We evaluate the production of large mass diphotons from quark annihilation at BNL 
RHIC and CERN LHC energies from central collisions of gold nuclei. The collision 
is assumed to lead to either a thermally and chemically equilibrated quark gluon 
plasma, or a free-streaming quark gluon gas having an identical initial entropy, or a 
chemically equilibrating quark gluon system, with the same entropy at T = Tc- We 
also obtain an estimate of hard photon pairs from initial state quark annihilation and 
find that the thermal production dominates the yield up to M ~ 4 GeV at RHIC, 
and up to 6 GeV at LHC. A simulation study of decay versus thermal diphotons is 
presented. 
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Relativistic heavy ion collisions are expected to lead to a confirmation of the 
QCD phase transition from hadronic matter to quark matter and an ephemeral 
formation of quark gluon plasma (QGP). Analysis of the S and Pb induced collisions 
at the CERN SPS have brought the search for QGP to an interesting stage. Many 
of the proposed signatures of QGP have been observed, though, most of them can 
also be explained to some extent in a purely hadronic picture. It is expected that 
the issue will be clearly resolved when we have results from BNL RHIC and CERN 
LHC, as much larger initial temperatures are likely to be attained there. 

Dileptons have long been considered as excellent probes of the early stages of 
relativistic heavy ion collisions. Large mass dileptons are likely to have their origin 
in the hot and dense stage of QGP and their detection would provide valuable infor- 
mation about this exotic state of matter. We suggest that an experimental detection 
of large mass diphotons can possibly provide a valuable confirmation of the results 
obtained from the measurement of dileptons. The theoretical understanding of the 
basic rates as well as the evolution of the interacting system has improved consider- 
ably since the early studies |||] in this direction, necessitating a re-evaluation of the 
importance of diphoton measurements. Thus, we have utilized the thermal mass of 
the annihilating quarks while evaluating the rates, which is more appropriate for 
quarks in a plasma. We have also considered evolution mechanisms which encom- 
pass the entire range of possible scenarios describing the expansion of the QGP. We 
give, as far as we know, the first estimate of hard QCD diphotons from the colliding 
nuclei. We also present a feasibility study of these measurements by performing a 
simulation of the likely measurements. 

Large mass diphotons from the QGP are produced from quark - antiquark 
annihilation. The cross section is given by |]^: 
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In the above A'c = 3, S = 1/2 and eq is the charge of the quark. This cross section 
diverges as rriq — > 0. The early works on diphotons mostly used mq = 5 MeV, 
which grossly overestimates the rates. We argue that as quarks in a heat bath 



acquire a thermal mass given by mth = VV^Tro^/S) T, this singularity can be 
regularized by making the natural choice, mq = mth- One can easily check that 
choosing m^ = 5 MeV, instead, will increase the basic cross-section by a factor of 
~ 3-5. Increasing mq from mth to 2mth, say, will decrease the cross-section by 
less than 50%. We shall see later that the explicit appearance of T in the basic 
cross-section here provides us with an additional probe for the evolution of the 
temperature |^. 

As remarked earlier, we investigate the production of thermal diphotons in 
three different scenarios, which essentially encompass the entire range of evolution 
dynamics one can imagine and thus give us the upper and the lower limits of the 
diphoton yield. As a first step, we consider the formation of a QGP in a thermo- 
dynamic (thermal as well as chemical) equilibrium at an initial time tq and initial 
temperature Tq in a central collision of two gold nuclei. This case has been studied 
in great detail in the literature. One can use the condition of isentropic expansion 
to relate the initial conditions to the particle multiplicity density (dN/dy) Q; 

dN 
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For a QGP consisting of u, d, and s quarks, and gluons, ag = 47.57r2/90. R is the ini- 
tial transverse dimension of the system. Assuming a rapid thermalization, one may 



get an upper limit on the temperature by taking tq = I/STq. Taking dN/dy = 1735 
for central collision of Au nuclei at RHIC we get an initial temperature of 478 MeV, 
while the corresponding numbers at LHC are 5624 and 860 MeV respectively |p. As 
we are interested in large mass photon pairs, we include only the emission from the 
QGP phase @, i.e., till the system reaches the phase transition temperature T^. We 
further neglect the transverse expansion which is known to affect the large invariant 
mass distribution only marginally Q. The proper time r^, when the temperature 
drops to Tc, is obtained from the Bjorken relation Tq = tq{Tq/Tc)'^ The invariant 
mass spectrum for diphotons will then be, 

dN tt/?^ I" '^1 
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At the other extreme, one may imagine that the system is formed at a proper 
time To with a temperature Tq after which the constituents free-stream away from 
the collision zone. One can evaluate the yield of diphotons from such a gas, in a 
manner similar to Ref. |5|, as 

dN _ .M^a{M,To)T^T^ 
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It is quite likely that the quark gluon system may be neither in thermal equilib- 
rium nor in chemical equilibrium when produced initially. Even though a thermal 
(kinetic) equilibrium may be attained quickly Q, the system may still be away from 
chemical equilibrium, which may evolve with passage of time due to a number of 



partonic reactions. Biro et al. [10|, for example, have studied the evolution of chem- 
ical equilibration with initial conditions determined from the HIJING model ||Tl[| . 
In order to have a meaningful comparison with the two descriptions above, we use 
the time variations of quark fugacity and temperature obtained by Strickland fl^ , 
with the stipulation that the particle multiplicity density at the end of the QGP 
phase is equal to the values of dN/dy used here. Now the diphoton spectrum is 
obtained from an expression similar to Eq. 4 with the additional introduction of the 



quark fugacity (Ag), taken from the Fig. 1 of Ref. |12]. 
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In Fig. 1 we give our results for the invariant mass distribution of the large 
mass photon pairs for RHIC energies. A number of observations are in order here. 
While the slope of the free-streaming description is essentially determined by Tq 
(see Eq. 5), the temperature drops from Tq to Tc in the hydrodynamic description, 
thus leading to an effectively steeper slope. Further, the yield for the equilibrating 
plasma is smaller due to the appearance of the square of the quark-fugacity in the 
expression and a more rapid cooling. Similar results are seen for the LHC energies 
as well (see Fig. 2). 

We must add here that the results for the equilibrated hydrodynamic evolution 
and those for the free-streaming are obtained by using tq = l/3To (see above), while 
those for the non-equilibrium hydrodynamics have tq = 0.3 fm/c. One may get yet 
another estimate for these measurements by starting the equilibrated hydrodynamic 
evolution from tq = 0.3 fm/c (see Figs. 1 & 2). Now the results for small M get 
quite similar to that for the non-equilibrium hydrodynamics. This looks surprising 
at first. Recall that our evolution scenarios are tailored to lead to a given final 
entropy (multiplicity). This ensures that during the later stages the conditions 
of the equilibrium and the non-equilibrium hydrodynamic evolution scenarios are 
nearly identical, leading to the similarity of these predictions at low M. 

The diphoton production from the hard QCD annihilation of quarks in the 
colliding nuclei at y = is given by; 



where x = M/y^ and ^/s is the cm. energy of the colliding nucleons. Np and Nt 
denote nucleons in the projectile and target respectively. We have introduced an 
arbitrary cut-off on the momentum transfer pc (=2 GeV), so that pQCD remains 



quark structure functions. Tab (= 293.4 fm ^) is the nuclear thickness function 
at zero impact parameter, corresponding to central Au-|-Au collisions. The result 
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applicable for such cases and used the MRSD— ' 



set of parametrizations |T^] for the 
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of this analysis for the RHIC and LHC energies are also given in Figs. 1 and 2, 
respectively. We find that thermal pairs dominate the spectrum up to 4 GeV at 
RHIC and up to 6 GeV at LHC. 

What additional information will large mass diphotons provide which we shall 
not already have from the measurement of dileptons? To explore this, we have 
evaluated the invariant mass distribution for dileptons as well as diphotons pro- 
duced from a thermally and chemically equilibrated QGP, which expands and cools 
to a temperature Tc (160 MeV), goes into a mixed phase where a phase transi- 
tion occurs, and when all of the QGP has converted into hadrons, again expands 
and cools till freeze-out temperature Tf (100 MeV) is reached. We consider a boost 
invariant longitudinal expansion without and with, a cylindrically symmetric, trans- 
verse expansion. We include the annihilation of pions in the hadronic matter and 
annihilation of quarks in the QGP as sources of dileptons [Q] and diphotons |]l|, ^. 

In Fig. 3 we have plotted the ratio of diphotons to dileptons as a function of M 
at RHIC and LHC energies. Firstly, we note that the structure seen in this ratio has 
its origin in the pion form-factor leading to the dileptons. The minimum at M ~ 0.8 
GeV comes from the p peak in the pion form-factor and the rapid rise there-after 
reflects the sharp decrease in this. The maximum around 1.4 GeV corresponds to 
the region where the (structure-less) quark annihilation process starts dominating. 
The differences in the 1-1-1 and the 3-1-1 dimensional cases are mostly confined to 
lower masses, which have their origin in the late stages of the evolution when the 
transverse expansion of the system is truly large. 

It is easy to verify that if we ignore the T dependence of the thermal mass of 
the quarks appearing in Eq. 1, a universal curve would be obtained for this ratio 
at large masses at all energies. The temperatures likely to be attained at LHC are 
larger than those at RHIC. Thus, in a purely longitudinal expansion these ratios at 
larger masses are smaller at LHC (recall the factor ~ ln(M^/2.7m^) in expression 
for cross-section given by Eq. 2). The transverse expansion of the system leads to 
a more rapid cooling of the plasma and thus higher temperatures are retained only 
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for a shorter duration. This will further increase these ratios, as indeed we see from 
Fig. 3. Thus we note that large mass diphotons sensitively reveal the details of 
the dynamics of the evolution of the plasma, especially when studied together with 
dileptons. The predicted variation, if observed, will also confirm our understanding 
of the temperature dependence of the thermal mass of quarks. 

Will it be possible at all to isolate diphotons from the huge background of 
photons coming from the decay of vr'^ mesons? It is quite clear that the background 
of decay photons can be arbitrarily reduced by choosing a large enough px- The 
search for diphotons will have to start with an accurate reconstruction of the decay 
photons so that any new source of photon pairs could be identified. We utilize 
the method of mixed-event analysis pioneered by the WA80 collaboration. We 
illustrate the situation at LHC energies and approximate the pT distribution of pions 
~ exp(— pr/T), with T ~ 300 MeV. The invariant mass spectrum for the photon- 
pairs from the decay of these pions in a sample of ten million central events is shown 
in Fig. 4(a). Similarly the invariant mass spectrum for the mixed events was also 
generated. This constitutes the back ground (see, e.g., Ref. [jl4|). A subtraction of 
the later from the invariant mass spectrum for the real events is shown in Fig. 4(b). 
The 7r° peak is seen very clearly. The diphotons from our calculations are seen to 
stand out clearly against this background. This simulation holds out the hope that 
diphotons could be isolated if high statistics data are available. It is important to 
note that the structure of the subtracted back-ground is quite different from the 
structure of the diphoton signal. Very similar results were obtained by fitting the 
background (e.g., in Fig. 4a) to a polynomial and subtracting the fitted spectrum 



from the invariant mass distribution |14]. 

In brief, we have evaluated the large mass invariant mass spectrum of thermal 
diphotons for a number of plausible scenarios in relativistic heavy ion collisions. An 
estimate of prompt diphotons is also made. It is shown that a comparison of dipho- 
tons and dileptons from heavy ion collisions can lead to an interesting confirmation 
of the temperature dependence of the thermal mass of quarks. Finally a simulation 
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of decay versus thermal diphoton production at LHC energies is performed which 
shows that such measurements could indeed be feasible, once high statistics data 
are available. 

We are grateful to Patrick Aurenche for useful communications in connection 
with the evaluation of prompt photon pairs. We are also grateful to Michael Strick- 



land for providing us tabulated values of T(t) and Xq{T) |12|. We thank Terry Awes 
and Itzhak Tserruya for useful suggestions in connection with the simulation study 
carried out here. 
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Figure Captions: 



Figure 1: Invariant mass distribution of photon pairs from qq annihilation for the 
Au+Au system at RHIC energies in a fully equilibrated QGP (Eq. Hydro), a free- 
streaming gas of quarks and gluons, and a chemically equilibrating (Noneq. Hydro) 
quark gluon system. Invariant mass distribution of hard photon pairs (Hard) from 
initial state qq annihilation for momentum transfer greater than 2 GeV is also shown. 

Figure 2: Same as Fig. 1 for LHC energies. 

Figure 3: Ratio of diphoton and dilcpton invariant mass distribution for a hadroniz- 
ing QGP at RHIC and LHC energies. Sources are qq 77 and qq IJ-'^IJ'~ in the 
QGP, and tt+tt" — 77 and tt'^tt" — /x"^//" in the hadronic phase. The solid and 
dashed lines indicate calculations with and without transverse expansion. 

Figure 4: (a) Distribution of invariant mass generated from photon pairs of the 
same events at LHC energies, (b) Invariant mass distribution of thermal diphotons 
and the mixed-event subtracted background of photon pairs. The difference is 
normalized so that the integral from M = 0.2 to 3.5 GeV vanishes. The hatched 
bar shows the 7r° peak. 
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